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ABSTRACT 


In May and June, 1983, NASA and USDA cooperated in the verification of 
a multi-sensor aircraft system developed to study soil moisture applications. 
This system consisted of a three beam push broom L band microwave radiometer, 
a thermal infrared scanner, a multispectral scanner, video and photographic 
cameras and an onboard navigational instrument. Ten flights were made of 
agricultural sites in Maryland and Delaware with little or no vegetation 
cover. Comparisons of aircraft and ground measurements showed that the system 
was reliable and consistent. Time series analysis of microwave and 
evaporation data showed a strong similarity that indicates a potential 
direction for future research. 
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1. INTRODUCTION 


A principal goal of the soil moisture remote sensing research is the 
evaluation and development of remote sensing technology for measuring and/or 
monitoring soil moisture. A secondary goal is to improve existing water 
management procedures through the use of this technology. In some cases, the 
adaptation of existing procedures is necessary because conventional soil 
moisture surveys are impractical and, therefore, the current water management 
procedures have developed under the assumption that this type of data could 
not be obtained. 

Truck and modeling experiments focus on the problems of defining optimal 
sensor systems and developing soil moisture estimation algorithms based on 
remotely sensed data. These studies allow evaluation of factors such as 
vegetation and soil type under well controlled conditions. 

Achieving the two primary goals of the project hinges on extrapolating 
these experiments to airborne systems, which are prototypes of future 
operational systems. Significant research questions, such as the effects of 
scene heterogeneity and instrument sensitivity at lower resolutions, can only 
be addressed using airborne sensors. 

Applications related objectives depend upon frequent large area coverage 
and rapid data turnaround. A dedicated airborne system was crucial to this 
line of research. 

Based on these considerations, NASA and USDA cooperated in the development 
of an optimal sensor aircraft system that would be available for extended- 
dedicated periods. It includes a multibeam L band microwave radiometer, a 
scanning thermal infrared system, and a mul tispectral scanner. All systems 
will be time and space registered and ultimately linked to an automatic 
geographic referencing system. 
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In 1983, this airciaft sensor system was assembled. The primary purpose 
of the study described in this report was to verify the performance of the 
system. Therefore, a fairly intensive ground sampling effort was included. 

A secondary objective was to evaluate the potential of repetitive observation, 

1 e., time series, in soil moisture and climate studies. 

Two agricultural test sites on the Delmarva Peninsula were selected for 
study due to their proximity to the research team and to the aircraft staging 
area. One check flight was conducted over the flight lines located in Maryland 
on March 14, 1983. Ten additional flights were made over the Maryland site 
from May 12 to June 27, 1983. Eight flights were conducted over the Delaware 
site during the same period. 

2. SITE DESCRIPTIONS 

A total of five flightlines were investigated, three near Carmichael, 
Maryland, and two in Kent County, Delaware. All of the sites are on the 
Delmarva Peninsula as shown in Figure 1. 

Maryland flightlines 1, 2, and 3 were located on or near the facilities of 
the Wye Research Center operated by the University of Maryland's Agricultural 
Experiment Station. Figure 2 is a black and white photograph of the area 
obtained on May 24, 1983, showing the field boundaries and flightlines. With 
the exception of March 14, 1983, all Maryland flightlines were flown west to 
east. Flightlines were flown at an altitude of 150 meters (500 feet) and line 

2 was also flown at 300 meters (1000 feet). 
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Figure 2. Maryland flight line 
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In March, all of the fields had corn or soybean stubble and/or weeds. 
During May and June a variety of bare soil conditions were present. Most of 
the fields were planted in corn, although, the tillage and mulch cover 
varied. Soils over this area were mostly silt loams. Additional details are 
presented in later sections. 

Flightline 4 was located near Sandtown, Delaware and was flown in hot* 
directions (east to west and west to east) at 500 and 1000 foot altitudes. 
Field boundarier* are shown in Figure 3. These fields included a variety of 
cover conditions over sandy loam soils. Figure 4 shows the location of 
flightline 5. Cover conditions at these test sites included corn and 
soybeans. Surface soil texture ranged from loamy sand to sandy loam. This 
line was flown in both directions at an altitude of 500 feet. 

3. GROUND TRUTH ACQUISITION AND PROCESSING 

3. 1 Soil Moisture 

Ground sampling at the Maryland sites was performed by personnel from USDA 
and NASA. Data at the Delaware sites were collected by personnel from 
Delaware State College. 

Gravimetric soil samples were collected using a 5 cm deep scoop with a 

3 

total volume of approximately 85 cm . The greatest number of samples were 
collected for the 0-5 cm soil layer. At the Maryland 'utes, a fewer number of 
samples were collected for both the 5-10 and 10-15 cm soil layer. For the 
Delaware sites the 5-10 cm soil layer was not sampled. 

All gravimetric samples were placed in metal cans and sealed with 
electrical tape. These were weighed at the end of the day to obtain wet 
weights. Samples were then placed in a drying oven at 104°C for 24 hours 
and reweighed to obtain dry weights. 
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Figure 3, Delaware flightline 4 and fields. 
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The number of samples collected in the Maryland fields depended on the 
field size; however, no fewer than six samples/f ield were obtained. Data were 
obtained on a grid centered on the flightline with a surface sample every 
150 m (500 feet). In most cases there were 12 points per field. Deeper depth 
samples were obtained near the field center at 2 or 3 locations. 

Several different schemes were used at the Delaware sites. These fields 
were sampled depending upon access and layout. All fields except number 17 
were sampled using a grid with a spacing of about 500 feet (150 m) between 
sanples. Due to access problems, only a border and one path in Field 17 were 
sampled. The number of samples per field averaged 16 for the 0 to 5 cm soil 
layer and 4 for the 10-15 cm layer. 

Gravimetric soil moisture values were averaged for each field on each 
date. This average was multiplied by the bulk density, described in the 
following section, to obtain the field average volumetric soil moisture values 
listed in Appendix I. 

An effort was made to collect the samples within 2 hours of the aircraft 
flights. In some cases, data from the irrigated Delaware sites had to be 
deleted because of irrigations that were being performed. 

On several sampling dates, especially at the Maryland sites, a significant 
portion of some fields was covered by standing water. Samplers were instructed 
that if a grid point happened to fall at a location with standing water, they 
were to note it as such in the field notebooks. The locations of the standing 
water followed very distinct patterns which in most cases coincided with local 
drainage. Comparing field average soil moistu ^ under these conditions to 
line samples of brightness temperature will require additional analyses beyond 
the detail presented in this report. 
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3.2 Bulk Density 

Determination of soil hulk density for the 0-5 cm soil layer was based on 
a volumetric displacement procedure that utilizes a specially designed bulk 
density ring with a hook gage and three one-foot-long bolts. The bulk density 
r 4 ng is placed on the ground and secured by driving in the three bolts. A 
sheet of plastic film is used to line the inside of the cylinder. Using a 
water-filled 500 ml graduated cylinder, the bulk density cylinder is filled to 
the hook gage and the quantity of water is recorded. This is returned to the 
graduated cylinder. The soil from the inside of the bulk density ring is then 
dug out to the desired depth of measurement and placed in a sealed container. 
The plastic liner is replaced and the ring is filled again. This amount of 
water is recorded. The soil wet and then dry weights are measured. The 
volume of soil removed is equal to the difference between the two water 
volumes used to fill the bulk density ring. The bulk density is computed by 
dividing the dry weight by this volume. 

Samples were collected several times over the study period at the Maryland 
sites, at least once in each field. Data were collected only once at each 
Delaware site after the completion of the experiment. Tour points were sampled 
in each field. 

All of the bulk density data collected over the experimental period were 
evaluated in terms of timeliness, tillage, cover, soils, and rainfall to 
estimate a bulk density value for every gravimetric sampling data. 
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Values assigned to each observation are included in Appendix I. In 

general, the following values were used: 

Condition Bulk density 

Maryland Delaware 



g/cm 

- 

Recently plowed 

1.11 

1.26 

Tilled, not recently plowed 

1.34 

1.36 

No till, idle 

1.48 

1.34 

Soil Texture 




Two sample cans were selected from one of the soil moisture sampling dates 
for each field. From each of these cans two samples were obtained (four total 
per field) to estimate the soil texture. 

The procedure used is based on a kit (Code 1067) produced by LaMotte 
Chemical Products Co. It yields the fractions of sand, silt and clay. Based 
on numerous repetitions using this procedure, we have found that it produces 
consistent results. 

Average texture results for each field based on the samples are listed in 
Table 1. Also listed in Table 1 are the USDA textural classifications based 
on sampling and the soils information derived from the soil surveys of the 
area (Soil Conservation Service, 1966 and 1971). With a few exceptions, the 
sampled textural classifications and the published soil textures were the same. 
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Volumetric moisture at specific tensions Sample bulk 
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3.4 Moisture-Tension Relationships 

The volumetric moisture content of the soils at specific tensions was 
determined using a pressure plate tensiometer. In such a device, a soil 
sample in an open ended ring is placed in a container to which pressure, in 
this case bottled nitrogen gas, is applied at a specific level. The water in 
the 9oil which is held by a tension greater than the applied pressure stays 
while the lower tension water is extracted through a porous plate. These 
samples are weighed before and after drying to determine the gravimetric 
moisture content. 

Samples were collected using a coring tool to obtain an undisturbed 
sample. Based on the volume of the ring, a bulk density value can be 
estimated. Moisture contents were determined for two samples per field at 
tensions of 1/3, 1, 3 and 15 bars. The 15 bar sample bulk densities must be 
estimated to obtain a volumetric moisture content because an unconsolidated 
sample must be used. 

The average values for each field are listed in Table 1. Figure 5 shows 
the results for three of the soils at the test sites plotted in the manner in 
which they are normally presented in the*literature. 

An alternative method for estimating moisture-tension values has been 
developed by Rawls et al . (1983). In this procedure, texture information, 
estimated from soils maps or sampling, is used with a series of graphs to 
estimate the volumetric moisture at specific tensions. These graphs were 
developed using an extensive data base. In this case, the sampled texture 
data were used with the graphs to estimate the 15 and 1/3 bar moisture 
contents listed in Table 1. In general, the results are similar to those 
obtained using the pressure plate tensiometers. The estimated values were 
higher, except for a few 1/3 bar samples. Considering the difficulties of 
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Figure 5. Laboratory moisture-tension relationships for 
selected Delmarva soils. 


pressure plate analyses, the problem of statistical sampling and spatial 
variations, the graphical procedure appears to yield consistent and reliable 
results. 

3.5 Soil and Surface Temperatures 

Soil temperature data were collected using metal-dial type temperature 
probes at depths of 5 and 15 cm. The 5 cm values were obtained at almost all 
soil moisture sampling points. Fifteen cm measurements were obtained at all 
15 cm soil moisture sampling locations, 3 or 4 near the field center line. 
Observed within field variations were typically less than 1°C at both depths. 

Surface soil/canopy temperatures were obtained using a hand-held thermal 
infrared radiometer. The sampling procedure involved the person walking from 
the edge to the center of the field and back and estimating the average 
response. All temperature data are listed in Appendix I. 

3.6 Cover Conditions 

The predominant crops on Delmarva peninsula are corn and soybeans. Corn 
is planted in early to mid May and soybeans in late May. Due to wet field 
conditions the soybean plantings were delayed this year. 

On March 14, 1983, all of the Maryland sites had stubble or weed cover. 
Figure 6a illustrates the density of corn stubble. Soybean stubble results in 
a higher percentage of cover, as shown in Figure 6b. 

In early May, most of the fields were planted in corn with the plants just 
beginning to emerge. By the beginning of June, the corn was between 10 and 
cm in height with less than 10X ground cover as shown in Figure 7. Heavy 
rains during this period made the soil surfaces fairly smooth in the tilled 
fields. At the end of June, some of the fields had close to 100% ground 
cover and some were just being planted with soybeans. 
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Several of the corn fields were net tilled, as a result, these usually had 
a ground weed cover and/or a heavy mulch or stubble cover. No till farming 
usually results in less soil surface evaporation due to the presence of the 
stubble. Figure 8 illustrates the ground cover in a corn field with a dense 
soybean stubble cover. Other fields had a weed cover for most of the study. 
Figure 9 illustrates the typical cover conditions. 

The Delaware flightlines included a wider variety of cover conditions. 
Fields 17 and 18 were vegetable crops and field 19 was mature barley. 

Observations during the study period were assembled and coded. For each 
field the following variables were estimated for each flight: 

1. Crop 

2. Percent ground cover 

3. Height 

4. Tillage 

5. Surface roughness 

6. Row direction 

Crop type codings that are used in Appendix I are listed in Table 2. 
Percent cover was assigned a value between 0 and 100%. Crop height was 
measured in centimeters. The tillage treatment codes are listed in Table 3 
and the visually estimated surface roughness codes are listed in Table 4. 

Row direction varied and for certain fields the reference row direction 
depends upon the flightline that was being considered. Codings used for the 
row directions are listed in Table 5. 


17 




ORIGINAL PAGE IS 
OF POOR QUALITY 




Table 2, n-op type codes. Tabl, Tillage condition codes. 


Code 

^rep i:ype 


Code Tillage condition 

i 

Bare soil 


< Tilled 

2 

Corn 


No till corn stubble 

3 

Soybeans 


i No till soybean stubble 

4 

Corn stubble 



5 

Soybean stubble 



6 

Mixed 



7 

Barley 



8 

Deciduous trees 



Table 

A. Surface roughness codes. 

Table 

5. Row direction codes. 

Code 

Roughness 

Code 

Direction relative to flightline 

i 

No recent tillage 

i 

Parallel 

2 

Smooth 

2 

Perpendicular 

3 

Moderate 

3 

Diagonal 

4 

Rough 

4 

None 


During the course of the study, a limited number of corn biomass samples 

were collected. On the average, 10 plants were removed, weighed, dried and 

weighed again to obtain the wet and dry biomass values per k lant. Vegetation 

water content was computed by subtracting the dry weight from the wet weight. 

The per plant values were then multiplied by the plant density of approximately 

6.6 plants per square meter to obtain the values on a square meter basis. 

Sample data are listed in Table 6. These values for an immature corn crop can 

be compared to data presented by O'Neill et al • (1983) for a mature canopy 

2 

which ranged from 5000 to 7000 gm/m for wet biomass. Based on vegetation 
effects models such as that presented by Jackson et al. (1982), the biomass 
observed on these fields should not have a significant effect on the 
sensitivity of the microwave brightness temperature to soil moisture 
variations . 
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Table 6, Corn biomass sample data. 


Field 

Date 

Plant 

height 

(cm) 


Per plant 


Per 

square meter 

Wet 

Biomass 

Dry 

Biomass 

Water 

content 

-g®- 

Wet 

Biomass 

Dry 

Biomass 

Water 

Content 

i 

6/15/83 

90 

132 

9 

123 

871 

59 

812 

i 

6 / 2 0 / 83 

120 

336 

27 

309 

2218 

178 

2040 

6 

6/i>/83 

30 

46 

4 

42 

304 

26 

278 

13 

6/15/63 

45 

76 

5 

71 

502 

33 

469 

13 

6/22/83 

60 

281 

28 

253 

1855 

185 

1670 

16 

6/15/83 

15 

13 

1 

12 

86 

7 

79 


Two of the fields that are reported in Appendix I were areas covered by 
deciduous trees. Field 15 was a cross section of the Choptank Fiver which was 
obscured by trees. Due to the heavy rains during the course of the study, the 
area covered by water varied. No ground verification of soil or water-covered 
area was obtained, only video data were used. Field 27 was a forested area 
that all three Maryland flight lines crossed. At times standing water was 
observed within this area. 

Field 17 was planted in peas in the early period of the study. The 15 cm 
rows made access difficult. A center pivot irrigation system wis present in 
this field. Field 18 was planted in several vegetable crops, including squash 
and cucumbers, and was partially bare. This field also had a center pivot 
system* The only other field with a sprinkler irrigation system was field 23* 
which was planted in sweet corn. 
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Field L9 had a mature barley cover over the entire period of study. It 
can be assumed that the water content of the crop was very low based on other 
studies involving winter wheat (O’Neill et al. 1983). Field 20 was 
discarded because the cover conditions were mixed corn and pasture. 

4. CLIMATOLOGICAL DATA 

Rainfall data were available at several weather station locations near the 
flight lines as shown in Figure 1. The Centerville, Denton, Royal Oak, Dover, 
and Milford data were extracted from the "Climatological Data for Maryland and 
Delaware" published monthly by NOAA. Delaware State College provided the 
rainfall data at the Darling and Willow Grove sites. 

The most complete and detailed record was from the class A station 
maintained by the Wye Agricultural Experiment Station (Wye Institute). This 
station included a weighing bucket raingage, evaporation pan, temperature, dew 
point, solar radiation, and wind speed. Data, excluding the raingage, are 
recorded in analog for:., and can be evaluated at 1-min intervals* The Wye 
Agricultural Experiment Station provided a 1-hour interval listing from which 
the daily values listed in Appendix II were computed for the period of May and 
June, 1983. Data for the Wye Institute raingage were obtained from strip 
charts. 

To illustrate the general meteorological conditions, Figure 10 was 
developed using the daily evaporation pan readings. This is a cumulative plot 
of evaporation minus rainfall, therefore positive values indicate more 
evaporation than rainfall. May started off with some very large rainfalls. 

The first few weeks of June were dry and at the end of June there were a few 
large storms. As described above, the daily evaporation pan readings reflect 
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both evaporation and rainfall. Hourly readings were reviewed to determine the 
separate rainfall and evaporation components. The estimated pan evaporation 
valves are listed in Appendix It. 

A detailed listing was not provided for the March 14, 1983, flight. 
ClimatoLogical records for the two weeks preceding the flight indicate that 
there was about 4 cm of rainfall with an event of approximately 1 cm on 
March 12, 1983. Daily air temperatures ranged from 35°F to 53°F on the day 
of the flight and the few days before. No freezing conditions were present 
within the preceding two weeks. 

Pan evaporation is directly related to the potential evaporation. A pan 
coefficient is used to relate the pan value to a lake evaporation estimate. 
This coefficient must be less than one and is usually about 0.7 (Jensen, 

1973). Also, it can be related to the potential evapotranspirat ion of 
particular crops. As described in Jensen (1973) the coefficient for an area 
like the Wye Agricultural experiment station, considering the immediate 
surroundings of the weather station, should range from 0.6 to 0.75. 

Potential evapotranspirat ion from grass can also be computed from the 
available meteorological data and the combination equation developed by Penman 
(Jensen, 1973), Using the average air and dew point temperatures, the wind 
speed and solar radiation listed in Appendix II, a daily potential 
evapotranspirat ion value was computed for each day. Results are listed in 
Appendix II. 

Pan evaporation and computed potential evapotranspiration were compared 
for the period of study. An optimization program was used to determine the 
value of a coefficient that would minimize the sum of squares of the residuals 
between the two. This analysis indicated that a pan coefficient of 0.65 was 
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best and had a coefficient of determination of approximately 0.92. These 
analyses indicate that the pan is performing as expected as an indicator of 
evaporation. 

5. AIRCRAFT AND SENSOR SYSTEMS 

The sensor system consists of an L band 3 beam push broom microwave 
radiometer, a thermal infrared scanner, a mul tispectral scanner, a 
photographic camera, a video camera and recorder, and a LORAN-C navigational 
system. 

5.1 Microwave Radiometer 

Additional details on the push broom microwave radiometer (PBMR) are 
presented in Table 7. It is a 3 beam L band system centered at 0° and 
+30°. The 3db (half power) points of the beams are located at roughly 
+10° of the beam center position. For the main beam this results in a swath 
(ground resolution) of 0.35 rimes the altitude. The two side beams will have 
swaths about 0.53 times the altitude. Beam centers for the side beams will be 
+0.58 times the altitude. Total coverage for the three beams will be a swath 
equal to 1.68 x altitude. At 150 m (500 feet) this will be 250 m (840 feet), 
at 300 m (1,000 feet) it is 500 m (1,680 feet), and 675 m (2,520 feet) at 450 
m (1,500 feet). Figure 11 illustrates the geometry of the PBMR system. 
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Table 7. L Band radiometer specifications. 


Antenna : 

Center frequency 

Bandwidth 

Polarization 

Beamwidth 

Beam efficiency 

Beam directions 

Radiometer: 

Type 

Center frequency 
Bandwidth 

Temperature sensitivity 
(short term stability) 

Accuracy (calibration stability) 
Operating temperature range 
Integration time (sec) 

Operation environ 


1,413.5 MHz 

30 MHz to 100 MHz, variable 
Horizontal 
20° at 3db points 
90 percent null to null 
0°, +30°, left and right of center 


Dicke 

1.413.5 MHz 

12.5 MHz 

T 0.7 °K, 1 sec integration time 
2°K 

+40°C to -40 °C 
0.5, sec 

unpressurized aircraft to 3,000 m 
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Figure U. Ground coverage and specifications for the push broo„ microwave radiometer 



5.2 Heat Capacity Mapping Radiometer (HCMR) 

The Heat Capacity Mapping Radiometer is a two channel scanning/ imaging 
radiometer. Channel 1 includes the spectral band from 0.55 - 1.1 microns and 
can be used to estimate albedo. The second channel provides a thermal infrared 
measurement over the interval of 10.5 - 12.5 microns. Figures 12 and 13 show 
the relative response of the two channels as a function of wavelength. 

The instantaneous field of view (IFOV) of the sensors can be up to 0.83 
milliradian (0.046°) ranging from + 45° off nadir. At an altitude of 300 
m (1000 feet) the nadir beam this would generate a ground resolution of 0.3 m 
(1 foot) resolution cell. A detailed description of the sensor system can be 
found in Bohse et al. (1979). 

The HCMR was operated at a scan rate that results in a gap of 
approximately 6 m between scans at an altitude of 300 m (1000 feet). At this 
rate, there would be approximately 50 scans per field which was considered 
adequate for this investigation. 

5.3 Ocean Color Scanner (OCS) 

The Ocean Color Scanner that was used in these studies is called UCS II 
and is a modification of an earlier system. OCS II has 12 channels as opposed 
to the 10 of the earlier version. 

Center wavelengths of the 12 channels are shown in Table 8. Unfortunately, 
spectral response curves for this system are not readily available. However, 
the curves for the 10 channel system are presented in Figure 14 to illus a^e 
typical conditions. 

The IFOV of the OCS can be 3.5 milliradians (0.2°) and it scans + 45° 
from nadir. At an altitude of 300 m (1000 feet) the footprint would be 1 m 
(3.5 feet). Additional details on this sensor system can be found in Blaine 


et al. (1977). 
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Processing 12 channels of high resolution data is impractical if rapid 
turnaround is desired. It is also unnecessary for this investigation. 
Therefore, only two channels will be processed for the biomass estimation 
procedures; channels 7 and 11. Ground coverage of this sensor is similar to t ; 
at of the HCMR. 

Table 8. - Center wavelengths of the OCS II. 



Center 

Channel 

wavelength 


i 

A25 

2 

464 

3 

509 

4 

549 

5 

587 

6 

631 

7 

669 

8 

711 

9 

756 

10 

793 

11 

828 

12 

863 


The photographic camera used was a 35 mm system with color film. Clock 
time was recorded on the image. 35mm color positive transparencies were 
generated for each flight. The primary image format used for ground-air time 
referencing were the color video tapes. These also had clock time on the 
image and were inexpensive. The aircraft on which the systems were mounted is 
a Short SC7 Skyvan, which is a twin engine plane. 


5.4 Data Processing and Calibration 

When the calibrated radiometer data were available for one-half second 
intervals, coverage for each field and run were extracted to compute an average 
and standard deviation for the brightness temperature. Start and stop times of 
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field boundaries were determined by reviewing the video tapes or photographs. 
At least one-half second of data was deleted from the beginning and end of the 
coverage to avoid minor variations related to aircraft parameters. 

In addition, coverage of each field was reviewed to determine if any mixed 
conditions resulted from flight line drift or changes in aircraft parameters. 
If these conditions were present, the data for that beam and field were not 
used . 


6.0 RESULTS 

Preliminary analyses of the data set focused on the relationship between 
the center beam (C) brightness temperature (T ) and the volumetric soil 
moisture (VSM) in the 0 to 5 cm soil layer. Each field was treated as an 
individual data set. 

Figure 15 illustrates the type of relationship that was observed for a 
tilled corn field (Field 1) with rows perpendicular to the flightline. The 
consistency of the response between different runs on the same day is very 
good. Measurements cover a wide range of soil moisture and brightness 
temperatures which adds reliability to empirical relationships. Some of the 
variation observed in Figure 15 can be explained by the increase in canopy. 

Data from day 173 had nearly 100Z corn cover and the brightness temperatures 
are higher than those we would expect from bare soils. 

Field 1 was one of the longer t .ae frames available and had a uniform set 
of cover and soil moisture conditions on each date. Therefore, this field 
will be used as a basis of comparison to other fields. 

Several other fields at the Maryland site were similar in soils and cover 
to field 1; these were 4, 8, 13, and 14. Graphical analysis of data from these 
fields showed that the patterns were very similar, although, the variability 
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was greater, as shown in Figure 16 for field 13. Regression parameters for 

brightness temperature and soil moisture listed in Table 9 indicate that these 

5 fields have similar intercepts and slopes. The coefficient of determination 
2 

R is quite a bit smaller in fields 8, 13, and 14, which reflects the vari- 
ability observed in the graphs. A major portion of the variability in these 
other fields can probably be attributed to differences between the average 
soil moisture of the field and the area actually covered by the sensor. 

Two of the fields, 2 and 7, remained in a bare soil condition through most 
of the experiment. However, they were tilled and planted for the last two 
flights. Regression results shown in Table 9 indicate similar parameters for 
the two fields. Results were also similar to those obtained for field 1. 

Three fields at the Maryland site, 3, 10, and 12, were not planted through 
most of the study but were covered with stubble and weeds. Qualitatively, 
field 10 had the densest cover, field 12 the smallest amount, and field 3 was 
in between. Plots of T. versus VSM revealed that the dynamic range of 
observed soil moisture vis related to the vegetation cover density, as 
illustrated by comparing Figure 17 for field 3 with Figure 16 for field 1. 

As expected from previous research, the observations of on these* fields 
tended to be greater than those observed over the corn field with less cover. 

The effect of the denser canopy of field 10 is apparent in Figure 18. 

These results are particularly interesting because on day P3 the field was 
plowed and we see thal the observation is very close to the field 1 regression. 
Once the canopy was removed, the field dries out faster and the data from the 
last flight show this to be true. 
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Figure 17. Brightness temperature vs soil moisture, 
field 3, center beam. 
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Figure 18. Brightness temperature vs soil moisture, 
field 10, center beam. 
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Table 9. Brightness temperature and 0-5 cm soil moisture linear regression results. 


TBC* = a+b(VSMl) 

Linear regression parameters 
TBL** = a+b(VSMl) 

TBR*** * a+b(VSMl) 

Field Intercept Slope 

Intercept Slope 

Intercept 

Slope 

(a) (b) 

(a) (b) 

(a) 

(b) 


i 


328 

-3.6 

.84 

334 

-4.2 

.88 

344 

-4.7 

.92 

2 


305 

-2.6 

.85 

293 

-2.5 

.89 

298 

- 2.7 

.81 

0 


259 

-1.1 

.44 

218 

.2 

.00 

219 

.2 

.02 

4 


221 

-3.4 

.76 

318 

-3.7 

.80 

262 

-1.8 

.35 

5 


317 

-2.6 

.32 

317 

-3.0 

.31 

321 

-3.0 

.32 

6 


345 

-4.3 

.69 

320 

-4.3 

.68 

318 

-4.2 

.61 

7 


319 

-2.9 

. 64 

303 

-2.8 

.61 

296 

-1.9 

.55 

S 


323 

-3.2 

.67 

315 

-3.4 

.68 

326 

-3.2 

.73 

9 


244 

-0.8 

.25 

255 

-1.4 

.44 

260 

-1.0 

.18 

10 


269 

-1.1 

.43 

248 

-0.7 

.18 

263 

-1.1 

.37 

11 


322 

-i.4 

.73 

321 

-3.6 

.78 

316 

-3.4 

.69 

12 


285 

-2.2 

.64 

283 

-2.5 

.66 

262 

-1.7 

.46 

13 


318 

-3.1 

.71 

325 

-3.4 

.78 

209 

-2.6 

.75 

14 

15 

16 


320 

-3.3 

.72 

322 

-3.7 

.85 

299 

-2.6 

.72 


227 

-0.6 

.11 

228 

-0.9 

.18 

204 

-0.0 

.00 

17 


343 

-6.4 

.75 

338 

-6.8 

.73 

343 

-6.4 

.77 

18 


334 

-4.8 

.85 

322 

-5.1 

.85 

335 

-5.2 

.83 

19 

20 
21 


261 

-1.4 

.75 

233 

-1.0 

.44 

265 

-1.8 

.18 


293 

-2.0 

.68 

284 

-1.8 

.57 

286 

-1.5 

.38 

22 


301 

-2.8 

.32 

296 

-3.2 

.51 

293 

-3.2 

.42 

23 


334 

-5.9 

.67 

322 

-5.8 

.65 

305 

-4.6 

.83 

24 


311 

-4.8 

.68 

309 

-6.3 

.91 

305 

-5.5 

.79 

25 


317 

-3.6 

.85 

299 

-1.9 

.62 

318 

-4.6 

.82 

26 

27 


308 

-2.6 

.87 

300 

-2.9 

.78 

300 

-2.4 

.78 

* 

TBC 

= Center 

beam brightness 

comparative 






** 

TBL 

= Left 

ii 

II 

n 






*** 

TBR 

= Right 

it 

II 

it 







Regression results from these fields are listed in Table 9, however, the 
limited range of observed soil moisture values makes them relatively useless 
on their own. 

One set of conditions that was observed in this study that had not been 
previously analyzed was mulched no till fields. Four fields, 5, 9, 16, and 22, 
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were no till corn with relatively dense ground covers of soybean mulch/ 
stubble. As a result of the mulch cover, there was less surface evaporation 
than on a tilled field; this resulted in a small range of observed soil 
moisture values. Figure 19 illustrates the typical relationship for field 5. 

A general interpretation of these results is that since the T_ exhibits a 
large dynamic range and the soil moisture does not, the soil moisture we are 
using is not indicative. Since there is a mulch layer present the moisture is 
fairly uniform with depth at the surface. This might suggest that the changes 
in T are the result of something related to the mulch. Perhaps there is a 

D 

build-up of moisture just under the mulch that affects the response. This 
aspect needs further study. 

Data collected over fields with lighter vegetation cover at the Delaware 
sites, such as field 23 shown in Figure 20, indicate a greater sensitivity of 
Tg to VSM. This result is typical of the change we expect between silt loams 
such as those at the Maryland sites, and the sandy loams at the Delaware sites. 

Field 19 had a mature barley cover through the series of flights. For all 
practical purposes, we can assume it to be the same cover on each day. Under 
these conditions, the effect of soil moisture changes under the canopy on the 
observed Tg can be evaluated. Figure 21 shows the plot of the data. There 
appears to be a distinct trend, although, the absolute level is not consistent 
with results from Field 1. 

Analysis of the data collected by the two side beams was performed through 
plots and regressions. Figure 22 shows the left and right beam measurements 
and the center beam regression for field 1. In this case there is no obvious 
difference between the two beams. Also, the regression from the center beam 
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data explains the side beam observations fairly well under dry conditions. 

Previous research, such as the study by Wang et al. (1983), have shown that 

there is little difference in the T_ between 0 and 30° incidence angles 

t> 

for bare soils under dry conditions, however, as the soil becomes wetter the 
difference increases. These results exhibit these same trends* 

Observations on the barley field (19) shown in Figure 23 suggest a uniform 
decrease in response for the side beams as compared to the center beam. These 
results must be interpreted with caution since there were no dry conditions 
observed. Comparisons to winter wheat data reported in 0*Neill et al . (1983) 
show that the expected decrease in T between 0° and 30° incidence 

D 

angles should be on the order of 15°K for this range of soil moisture 
conditions, which is approximately the amount observed here. Regression 
analyses of the side and center beams, which are summarized in Table 9, show 
that there is a strong linear relationship on a field basis with the exception 
of the no till and weed/stubble covered field 5. In addition, in most of 
these same fields there were poor linear relationships for the left and right 
beams. These results once again suggest spatial soil moisture patterns in 
these fields. 

Due to the frequency of the flights in Hay and June and the variation in 
weather conditions, the data collected in this study are suitable for the 
analysis of temporal patterns of T and their relationship to climatic 

D 

variables beyond surface soil moisture. For the present time, the analysis 
will use the cumulative evaporation pan values collected at the Uye 
Agricultural Experiment Station. This value represents an approximation of 
the cumulative outflow and inflow to the soil column. Pan evaporation is 
usually greater than free water or lake surface evaporation which is used to 
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Figure 23. Brightness temperature vs soil moisture, 
field 19, side beams. 
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approximate the potential evapotranspiration. Potential evapotranspirat ion is 
used in conjunction with soil and vegetation parameters to compute actual 
evapotranspiration. 

The cumulative daily evaporation pan data (EE^) will be positive when 
evaporation for the period exceeds rainfall and vice versa. It was set to 
zero on Julian day 124 and plotted in Figure 10. 

Figure 24 includes a plot of T as a function of Julian day. On its own 

D 

it is difficult to interpret, however, when it is plotted with the ^E^ 
relationship using an arbitrary scale, a remarkable similarity of pattern is 
apparent. The T n is responding in a nearly identical way as the F E . During 

B p 

the extended drying period in early June the relationship is very strong. This 
might be expected with a bare soil condition. However, similar patterns are 
observed between the last two flights when there was nearly 100% cover. 

With these results in mind, data from several other fields were analyzed. 
Figure 25 shows the results obtained for field 3 (tall weeds). Here the 
dynamic range is ouch less, however, the pattern is still apparent although 
not as dramatic. 

Figure 26 is a plot of the data collected over field 5. This response is 
between the two extremes of fields 1 and 3, which is expected based on its 
cover. What is very important here is the fact that the T^ values follow a 
very clear pattern which, based on Figure 17, doesn't necessarily correspond 
to the VSM in the 0 to 5 cm soil layer. 

Data collected over the two forest covered fields, 15 and 27, did not 
include ground conditions. Therefore, the only comparison that can be made is 
to general climatic variables. Having shown that there is reason to believe 
that a relationship exists between the Tg and I F^, comparisons on the 
forest sites is possible. 
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Figure 27 is a plot of the data collected over field 27. A general trend 
similar to those observed in the weed field is present, however, the day-to-day 
variation is small and on the same order as the variation on a given day. 

The results over field 15, the Choptank River cross-section, are quite 
different. As shown in Figure 28, the pattern is quite distinct. It is 
hypothesized that the reason for the variation in T is the extent of the 
flooded area beneath the trees. As the streamflow receded the width of a 
shallow channel such as this section of the Choptank River also recedes. 

Another aspect of the forest sites that was considered was what difference 
might be expected in between beams with such a dense cover. Responses for 
field 27 showed no consistent difference on any day. Greater variability was 
observed for field 15, which is probably related to spatial patterns. 

7 . SUMMARY 

In May and June, 1983, NASA and USDA cooperated in the verification of a 
multi-sensor aircraft system developed to study soil moisture applications. 

This system consisted of a three beam push broom L band microwave radiometer, 
a thermal infrared scanned, a raul tispectral scanner, video and photographic 
cameras and an onboard navigational instrument. Ten flights were made of 
agricultural sites in Maryland and Delaware with little or no vegetation cover. 
Comparisons of aircraft and ground measuremnts showed that the system was 
reliable and consistent. Time series analysis of microwave and evaporation 
data showed a strong similarity that indicates a potential direction for 
future research. 
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Brightness temperature vs Julian day, field 27, center beam. 


5 





8. ACKNOWLEDGEMENTS 


The authors wish to express their sincere appreciation to the personnel of 
the Hydrology Laboratory, USDA/ARS/BARC , the Hydrological Sciences Branch 
NASA/GSFC and the Department of Agricultural and Natural Resources of 
Delaware State College for their assistance in the data collection effort* We 
are especially thankful to Dr. Russ Brinsfield of the Wye Agricultural 
Experiment Station for his assistance and cooperation in this project and to 
the individual farmers who allowed us to sample their fields. 


9. REFERENCES 

Blaine, L. R. , C. J. Murphy, III, and W. L. Barnes, 1977. The Ocean Color 
Scanner Experiment, GSFC X-Docuiment X-941-77 — 153, NASA, Goddard Space 
Flight Center. 

Bohse, J. R., M. Bewtra, and W. L. Barnes, 1979. Heat Capacity Mapping 
Radiometer (HCMR) Data Processing Algorithim, Calibration, and Flight 
Performance Evaluation, Tech. Memo. 80258, NASA, Goddard Space Flight 
Center. 

Jackson, T. J., T. J. Schmugge, and J. R. Wang, 1982. Passive Microwave 

Remote Sensing of Soil Moisture Under Vegetation Canopies, Water Resources 
Research, Vol. 18, pp. 1137-1142. 

Jensen, M. E. (ed), 1973. Consumptive Use of Water and Irrigation Water 
Requirements, ASCE, New York, NY. 


53 



O'Neill, P., T. Jackson, B. Blanchard, R. van den Hoek, W. Gould, J. Wang, 

W. Glazar, and J, McMurtrey, III, 1983. The Effects of Vegetation and 
Soil Hydraulic Properties on Passive Microwave Sensing o± Soil Moisture: 
Data Report for the 1982 Field Experiments, NASA TM 85106. 

Rawls, W. J. and D. L. Brakensiek, 1983. A Procedure to Predict Green and 
Ampt Infiltration Parameters, ASAE Winter Meeting, Chicago, Illinois 
Soil Conservation Service, 1966. Soil Survey Queen Annes County, Maryland, 
U.S. Government Printing Cfiice, Washington, D.C. 

Soil Conservation Service, 1971. Soil Survey Kent County, Delaware, 

U.S. Government Printing Office, Washington, D.C. 

Wang, J. R. , P. E. O'Neill, T. J. Jackson, and E. T. Engman, 1983. 

Multifrequency Measurements of the Effects of Soil Moisture, Soil Texture, 
and Surface Roughness, IEEE Trans, on Geoscience and Remote Sensing, 

Vol. GE-21 , No. 1, pp. 44-51. 


54 



APPENDIX I 

AIRCRAFT AND GROUND MEASUREMENTS 
DELMARVA SITES, 1983 


LISTING CODES 

Date - Calendar date 

Day = Julian day, 1983 

Field » Field number 

Line = Line number 

Run = Run number 

Dire = Flightline general direction; 0 = west to east and 1 = east to west 

Crop = Crop type code, se^ able 2 

PCover - Percent ground cover 

Hgt = Crop height (cm) 

Row = Row direction, see Table 5 

Till = Tillage code, see Table 3 

Rough = Soil surface roughness code, see Table 4 

BD = Soil bulk density (g/cm^) 

Tl = Soil temperature at 5 cm (°C) 

T2 = Soil temperature at 10 cm (°C) 

T3 = Soil temperature at 15 cm (°C) 

TIR = Field thermal infrared temperature, hand held (°C) 

VSM1 = Volumetric soil moisture 0-5 cm (%) 

VSM2 = Volumetric soil moisture 5-10 cm (%) 

VSM3 = Volumetric soil moisture 10-15 cm (%) 

SD1 = Standard de'* iation of volumetric soil moisture 0-5 cm (%) 

SD2 = Standard deviation of volumetric soil moisture 5-10 cm (%) 

SD3 = Standard deviation of volumetric soil moisture 10-15 cm (%) 

TBL = L band brightness temperature, left beam (°K) 

TBC = L band brightness temperature, center beam (°K) 

TBR = L band brightness temperature, right beam (°K) 

SDL = Standard deviation of brightness temperatures, left beam (°K) 

SDC = Standard deviation of brightness temperatures, center beam (°K) 

SDR 55 Standard de /iation of brightness temperatures, right beam (°K) 
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APPENDIX II 

CLIMATOLOGICAL DATA - DELMARVA PENINSULA 
May and June 1983 


LISTING CODES 


DAY = Julian Day, 1983 

R1 = Daily rainfall, Wye Institute, MD (cm) 

R2 = Daily rainfall, Centerville, MD (cm) 

R3 = Daily rainfall, Denton, MD (cm) 

R4 = Daily rainfall, Royal Oak, MD (cm) 

R5 = Daily rainfall, Dover, DE (cm) 

R6 = Daily rainfall, Milford, DE (cm) 

R7 3 Daily rainfall, Willow Grove, DE (cm) 

R8 = Daily rainfall, Darling, DE (cm) 

EVAP = Daily pan evaporation (cm) 

TAIR 3 Mean daily air temperature (°C) 

TDEW 3 Mean daily dew point temperature (°C) 

WIND = Daily wind (km/day) 

RAD 3 Solar radiation (Langleys) 

PET = Penman daily potential evapotranspiration (cm) 


PRECEDING PAGE BLANK NOT FILMED* 
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